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Traveling waves in two-component Bose-Einstein condensates whose dynamics is described by the 
Manakov limit of the Gross-Pitaevskii equations are considered in general situation with relative 
motion of the components when their chemical potentials are not equal to each other. It is shown 
that in this case the solution is reduced to the form known in the theory of motion of S. Kowalevski 
top. Typical situations are illustrated by the particular cases when the general solution can be 
represented in terms of elliptic functions and their limits. Depending on the parameters of the 
wave, both density waves (with in-phase motions of the components) and polarization waves (with 
counter-phase their motions) are considered. 
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I. INTRODUCTION 

Realization of Bose-Einstein condensates (BECs) of 
atoms which can occupy several quantum states at ex¬ 
tremely low temperatures has drawn interest to nonlin¬ 
ear dynamics of such multi-component systems (see, e.g., 
the review article Q)- In particular, if the components 
can mix, then in such a condensate two types of linear 
waves can propagate—usual sound waves (that is, den¬ 
sity waves) with in-phase oscillations of the components 
and so-called polarization waves with counter-phase os¬ 
cillations. Correspondingly, generally speaking, there ex¬ 
ist two Mach cones and two channels of Cherenkov ra¬ 
diation what leads to considerable changes in the char¬ 
acter of excitations in the system compared with one- 
component situation. The same holds true for the non¬ 
linear excitations—solitons and breathers. For example, 
oblique solitons can be generated by the flow of two- 
component BEC past a non-polarized obstacle which re¬ 
pels both components Q and oblique breathers are gen¬ 
erated in the case of polarized obstacles which repel one 
component and attract the other one Q. 

The dynamics of two-component BECs becomes even 
richer if one component moves with respect to another. 
As was found experimentally in [3^^ j the relative motion 
of components leads to generation of nonlinear periodic 
waves of polarization. In these experiments the com¬ 
ponents correspond to the quantum states |I,—I) and 
12,—2) of the hyperfine structure of ®^Rb atoms with 
very close values of their inter-atomic scattering lengths. 
Hence their quasi-one-dimensional dynamics in elongated 
cigar-shaped traps can be described with high accuracy 
W the Gross-Pitaevskii equations of the Manakov type 
[7|. In the standard non-dimensional units these equa¬ 
tions can be written in the form 

*V'fe,t + ^V'fc,a:a:-(|'01P+|'02n'0fc =0, k = l,2, (1) 

where 'ipk denotes the wave function of the fcth compo¬ 


nent, X is the coordinate along the trap, and t is the 
time variable. Such multi-component (vector) nonlinear 
Schrodinger equations appeared first in nonlinear optics 
Q and they have been studied intensely in this physi¬ 
cal context, where it is natural to suppose that the wave 
numbers (which are analogues of velocities of the BEC 
components) of both components are equal to each other 
(see, e.g., Q), but the interaction constants are differ¬ 
ent. Thus, so far the situation with equal nonlinearity 
constants and non-vanishing relative velocity of the com¬ 
ponents was studied very little. An important particular 
case of such a type of solutions is so-called dark-bright 
soliton with vanishing of one of the background densi¬ 
ties far enough from the soliton location. This means 
that the component with non-vanishing background den¬ 
sity forms a “trap” for another component localized in¬ 
side such a trap (see, e.g., [a El)- Although this so¬ 
lution of the Manakov system describes an important 
type of nonlinear excitations in two-component BECs, it 
cannot explain dense lattices of dark-bright solitons ob¬ 
served in recent experiments An attempt of such 

an explanation was done in Refs. [l2l - [l5l| where partic¬ 
ular cases of nonlinear waves with relative motions of 
the components were studied. Indeed, solutions in the 
form of counter-phase oscillations (polarization waves) 
were found in these papers, however, they were limited to 
BECs with equal chemical potentials and this condition is 
quite restrictive for adequate description of experimental 
observations. 

In this paper we shall consider the general situation 
with non-equal chemical potentials for the case of one- 
phase traveling waves. It will be shown that in this case 
the Manakov system can be reduced to the equations 
studied first by Sophie Kowalevski in her theory of ro¬ 
tation of the so-called Kowalewski top [T^ . [T^ . (Similar 
reduction was performed for the Manakov system with 
attractive (focusing) nonlinear interaction and without 
relative motion of the components in Refs. [Tsl - f^ .) 

The physical conditions that the densities of the com- 
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ponents must be positive and non-singular impose heavy 
restrictions on the admissible solutions of the Kowalevski 
equations. The typical situation will be illustrated by 
several particular cases. In particular, it will be shown 
that the solutions studied previously in [IMl can be 
obtained as a special limiting case of the general so¬ 
lution of the Kowalevski equations. Another particu¬ 
lar case of the dark-bright solitons is also obtained as 
a result of simple degeneration of the Kowalevski equa¬ 
tions. The so-called Appelrot-Delone class of solutions 
of the Kowalevski equations, studied previously in the 
context of rotations of the Kowalevski top, leads now in 
the context of two-component condensate flows to the 
nonlinear density waves with very specific dispersion law 
conditioned by the relative motion of the components. 
In the soliton limit this periodic solution reduces to the 
dark-dark solitons which generalize the known Manakov 
soliton solution to the situation with relative motion of 
the components. The general nonlinear wave in the two- 
component BEC is illustrated by the so-called Legendre- 
Jacobi case when the hyperelliptic integrals are reduced 
to the elliptic ones. The physical implications of the 
found solutions are discussed in the conclusion of the pa¬ 
per. 


II. EQUATIONS OF MOTION, THEIR 
INTEGRALS AND GENERAL SOLUTION 


where U is a constant velocity of the wave. Let us in¬ 
troduce the imaginary “time” variable t = —2i^. Then 
after obvious integrations we get 

dpk dwk 1 2 , , o 

-^ = ak- PkWk, -^ =+ Pi + p2 - Pk, ( 6 ) 

where ak, Pk are the integration constants and Wk = 
Vk - V. 

It is worth noticing that integration of the first pair 
of Eqs. ([3]) under ansatz (O gives expressions for Uk in 
terms of pk, 


Ufc(a = ^+^, * = 1,2. (7) 

Similar integration of the second pair of Eqs. ([3]) followed 
by exclusion of Uk with the use of Eqs. 0 yields 







Pk- 


( 8 ) 


As follows from these relations, the constants ak, Pk, 
k = 1, 2, are real. In case of a uniform flow with con¬ 
stant Pk = Pko and Uk = Uko we have Pk = pio + P 20 + 
'^fc/(2Pfco)- Haace, the parameters Pk are related with 
the chemical potentials pk = ul/2 -|- pio -I- P 20 by the 
formulae 


It is convenient to transform the Manakov system ([T]) 
to the hydrodynamic-like form by means of the Madelung 
transformation 


ipk = y/pk exp 




( 2 ) 


where pi ,2 are constants. In a standard way we arrive at 
the system 


Pk,t T {PkUk^x — 0, 


'^k,t 



Pk,a 


'ipk 


= 0 


(3) 


1 o Vak „ , 

pk = -zV H-f Pk, k = 1,2. (9) 

z PkO 

The solution studied in Refs. (l2l - [T5l| were limited to the 
case Pi = P 2 what meant that in the uniform case the 
chemical potentials are equal to each other: pi = p 2 . 
Here we will consider the general case including situations 
with Pi ^ P 2 - 

Our approach is based on the fact that the system ([6|) 
is Hamiltonian with the Hamiltonian 

H + P 2 W 2 ) - ^{pi + P2)^ 

-I- aiwi -I- a2W2 + Pipi + P2P2 


with real variables. Here pi ^2 denote the densities of the 
components, ui ^2 denote their flow velocities, and pi ^2 
are their chemical potentials. 

The complex substitution Vk = Uk — ipk,x/ {2-pk) casts 
the system ([3]) to a mathematically simpler form 


Pk,t + ( PkVk + -^Pk,x j — 0, 

Vkp + ( + Pi + P 2 — 2^k,x ) = 0. 


(4) 


A traveling wave is described by one-phase solution of 
Eqs. dH) with all the variables depending on ^ = x — Vt 
only. 


Pk=Pk{0: Vk=Vk{0, ^ = x-Vt, (5) 


and Poisson brackets 

{pi,Pj} = {w^,Wj} = 0 , {wi,pj} = S^J. 

The corresponding equations of motion 

dpk OH dwk dH 

dr dwk ’ dr dpk 

possess the integral of energy 

H{pk,Wk) = h = const. 


( 11 ) 


( 12 ) 


(13) 


For complete integrability of this system with two degrees 
of freedom we need, according to the Liouville-Arnold 
theorem (see, e.g., [21|), one more integral. One can 
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check that there is such an integral quadratic in momenta 

Wk- 


K = — pip 2 {wi — W 2 f + 2(wi — W 2 ){aip 2 — a 2 Pi) 

- (/3i - I32)[piwl - P 2 wl + pI- pI 

- 2{aiwi - 02 ^ 2 ) - 2(/3i/9i - /32P2)]- 

(14) 

Thus, integration of the system dm) can be reduced to 
quadratures. 

If /3i ^ P 2 , then we can make a canonical transforma¬ 
tion 


_ (gi + /3)(g2 + /?) _ (gi - /3)(g2 - /3) 

2/3 ’ 2/3 

{qi - P)pi - {q 2 - I3)p2 

Wi = - 

qi - <?2 

(li + P)Pi - (12 + I3)p2 

W2 = -, 

qi - <?2 


(15) 


where we have denoted /3 = /3i — /32 (the limit /3 —>■ 0 
will be considered below). As we shall see, the dynamics 
is separable in these new variables (/i, Pi, i = 1, 2. The 
Poisson brackets preserve their canonical form 


= {PoPj} = 0, {Pi,qj} = (16) 


and the Hamiltonian becomes 
„ ^ ( 9 ? - P‘^)pI - 2[(ai + a2)qi - (ai - a2)l3]pi 
2(92 - 91 ) 

(gj - P‘^)pI - 2[(ai + <^2)92 - jai - a 2 )l 3 ]p 2 ny) 
2(91 - 92) 

- +gi + 5192 - (/3i +/ 32 )(gi -f 52 ) -/3^]. 


The equations of motions are given by 
dqi dH 
dr dpi 

_ {qf - I 3 ‘^)pi + {ai - Q;2)/3 - (qi + 02)51 

gi - 92 

dpi dH 

dr dqi 

_ (pi -P2)[iai - Q2)/3 - (ai + 02)52] 

(51 - 52)^ 

_ (9i + d'^)pl - (gj - I 3 ‘^)pI + 29i92pj 

2(gi - 52)^ 

+ + d2- 251 - 92), 


(18) 


and similar equations can be written for 92 and 
P 2 - They have two integrals of motion—the energy 
H{qi,pi,q 2 ,p 2 ) = h = const and 


K{qi,Pi,q2,P2) = {2(ai -f a2)(Pi -P2)5i52 

- (P?gi -P252)5ig2 + /3 ^(p?52 -P25i) 

- 2/3(ai - a 2 ){piq 2 -P 2 gi)}/( 5 i - 52 ) 

+ (/3i + /32 - 51 - 52)5152 = k = const. 


As was mentioned above, integration of Hamiltonian 
system with two degrees of freedom and two integrals of 
motion can be reduced to quadratures. Actual integra¬ 
tion can be performed in our case as follows. Eliminat¬ 
ing variables 92 , P 2 from the integrals H{q,p) = h and 
K{< 1 jP) = k, we obtain 

$ = ( 5 ?-^^)(^i+/32-51-P?) , 20 ) 

-I- 2pi[(ai -I- 02)51 - ,5(01 - ^ 2 )] - 2/191 + k = 0, 


what demonstrates the mentioned above separation of 
variables. Taking into account Eqs. ra, we get 


5$ 

dpi 


2(9?-/32)(/3i+/32-9i-P?) 


2[/3(oi - 02 ) - (oi -I- 02 ) 91 ] 


-2(91 



( 21 ) 


We solve equation (EH)) with respect to pi and substitute 
the result into (1211) . After this and similar manipulations 
with the variables 92 , P 2 we obtain the system 

±i/-TZ{qi) = -(91 - 52 )^, 

,_ ( 22 ) 

±i/-TZ{q2) = (91 - 52)^, 
dr 

where 


n{q) =9^-(/3i+/32)9"-2(/32-h)93 

— [(^^i + < 112 )^ — 2(/3i -l- /32)/3^ -l- k\q^ 
+15[/5^ — 2/i/3 -l- 2(oj — 02)]9 

— i5^[(/3i -|- ^ 2 ) 0 ^ — k + (oi — 02 )^] 


is a 5th degree polynomial with respect to 9 . Then after 
simple manipulations we arrive at the system 

dqi dq 2 


qidqi 


52 1/52 


\/^(5i) ^"^( 52 ) 


= 0 , 


= ±2d^. 


(24) 


where we have returned to the real variable ^ = jr/ 2 . 
Sometimes it is convenient to rewrite this system in the 
Kowalevski form 


dqi _ 20 n{qi) dq 2 _ 20 H{q 2 ) 

di qi - q 2 qi - q 2 

The systems (Eli) or (E51) can be solved formally in 
terms of Riemann 0-functions (more precisely, in terms 
of Gopel and Rosenhein hyperelliptic functions; modern 
exposition of this method can be found, e.g., in ( 2 ^) but 
such a form of the general solution is mathematically in¬ 
volved and hardly can produce essential understanding 
of physical behavior of waves in a two-component BEC. 
Therefore we shall confine ourselves here to the most im¬ 
portant particular solutions which provide useful infor¬ 
mation about such typical nonlinear excitations in BEC 
as periodic waves and solitons. 

























III. NONLINEAR WAVES IN A 
TWO-COMPONENT EEC 


A. Limit /3 —^ 0 
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The physical variables pi and p 2 (i-e. densities of EEC 
components) must be positive and this condition imposes 
important restrictions on the variables qi and q 2 which 
obey the systems (|24l) or (l25ll . Supposing for definiteness 
that /3 > 0 , it is easy to find that qi and q 2 can vary in 
the intervals (see Fig. 1) 


If ^ = /3i—/32 —>■ 0, then we must have i^i ^ 0 and ^2 —>■ 
0 to satisfy the conditions (l26l) . At the same time, to get 
in this singular limit from Eqs. (HB finite values of pi 
and p 2 , we have to define new variables and parameters 
as 

qi = /3qi, vi = I3vi, 1^2 = 13^2 (30) 


-13 <qi< 13, q 2 > (3, 
or qi >13, -P <q 2 < /3. 




(26) 


so that Eqs. m reduce to 

Pi = ^ 92 ( 1 + 9i), P2 = ^< 12(1 - Qi)- 


(31) 


The Kowalevski equations (1^51) are then transformed to 


dqi 2 ,^ 1 / 31 / 41/5 /— -——-— 

-JT = -V (91 “ - 92), 

92 


(32) 


” 1 

p 

---1 

_ _ _ _ _ _ 


-p 


p 


91 


-P 


FIG. 1: Regions of variations of the parameters q\ and 52 (see 
Eq. (1261) 1 corresponding to the conditions of positivity of the 
densities pi and p 2 defined by Eqs. jm. 


Formal integration of the system 
n dq Z-® dq 


' qio 
ri 


^/W) 

qdq 


+ 


+ 


' <120 
rQ2 


'920 


qdq 

^/W) 


yields 

= 0 , 

= ±2e, 


(27) 


where gip and 520 are integration constants equal to the 
values of qi and 52 at ^ = 0, respectively. Every solution 
of the system (IMl) is parameterized by five zeroes Vi, 
f = 1,..., 5, of the polynomial 
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nq) = 11(9 - 1'^) (28) 

i-l 


We notice that the expression p = pi + P 2 = 9i + 92 
reduces to p = <72 in this limit. Introducing also gq = 
cosP, i7i = cos 6 * 1 , V 2 = cos 02 , = pi, Vi = P 2 , = P3, 

we arrive at the equations 


dcosO 

dp 

d^ 


- y/PlP2P2 ^_ cos0i)(cOS02 — COS0), 

P 

-2\/(P - Pl)(P2 - P)(P3 - P) 

(33) 


identical to the equation obtained for this special case in 
Refs. [iMI. 

It is worth noticing that the integrals of motion (fT^ 
and m can be cast after introduction of these new vari¬ 
ables p and 9 to the form 




8 p 


El 

2 


■ Pp 


p^ sin^ 0 • 0| -I- 8 [ai -I- ~ “ ^ 2 ) cos0] 

8 p sin^ 0 


= h. 


(34) 


p^ sin^ 0 • 0| -I- 8 [af -I- — {a\ — a\) cos 0] 

4ih?0 (35) 

— (ai -I- 02 )^ = k. 


and, depending on their values, we obtain different 
classes of solutions. Since the solution is symmetric with 
respect to transposition of gi and 52 , for definiteness we 
assume that they change in the intervals (the zeroes Vi 
are numerated in the order of their values) 

vi < qi < V 2 and 1/3 < 92 < (29) 

where the polynomial 72.(g) is positive. 


The angle 0 can be excluded from Eq. (IMl) with the use 
of Eq. (l35l) what gives the equation for a single variable 
P, 

p| = 4p^ - 8/3p^ -I- 8hp - A[k + (oi -I- 02 )^] (36) 

which is another form of the second equation (1551) . Its 
solution can be expressed in terms of elliptic functions. 
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When p is known, then 9 = 9{^) can be obtained by 
integration of the equation (1551) or 



k + {ai + a 2 )^ 
2 


af + a 2 — (af — a^) cos 6 
sin^ 9 


(37) 

which also coincides up to the notation with the first 
equation (l351l . When their solutions are found then the 
components densities are given by (1311) transformed to 


Pi(0 = PiO cos^ P2{0 = piO sin^ ^ (38) 


More details about these solutions can be found in 
Refs. (T2l - [l5l| . In particular, they include also the well- 
known Manakov dark-dark soliton solution. 


B. Appelrot-Delone class of solutions 


As we shall see, the second choice cannot give the soliton 
solution, so we shall consider the first one. 

In standard way we obtain 

qiiO = i^i + (;^2-J^i)sn2(\/I^^7717;'(^-^o),m), (43) 

where 


m = 


V2 - vi 

V 3 -V 1 


(44) 


and, to simplify the notation, from now on we shall put 
the integration constant ^0 = 0- Then the components 
densities are given by 


Pi 

P2 


[P + I'l + {i '2 - I'l )sn 2 WI's - I'l m)], 
[P - 1^1 - {1^2 - i 2 i )sn 2 (V I's - 121 C, m)], 


(45) 


The systems (l24)) and (|25l) were applied for the first 
time to a real mechanical problem by Sophie Kowalevski 
in her theory of rotation of the so-called Kowalevski top 
[ 1 ^ . After that some particular especially remarkable 
motions of this top were discussed by other authors, in 
particular, by G. G. Appelrot and N. B. Delone (see, e.g., 
(T^l. Here we shall apply their method to the special case 
of nonlinear motion of a two-component BEG which we 
shall also call the Appelrot-Delone case. 

Let us suppose that the polynomial 'R-{q) has a double 
root q = V, the other roots we denote vi < V 2 < vs, that 
is we have 


7 ^( 9 ) = {q-vf{q-vi){q-V2){q-V3) = {q - vfni{q), 

(39) 

where TZi {q) is the 3rd degree polynomial with the roots 
1 ^ 1 , V 2 , t's. In this case it is convenient to use the system 
(1251) written now in the form 


2 (gi - v)\/ni{qi) = {qi - 92)^, 

dS. 

2(92 - l^)\/^u(^ = -(91 - 92 )^- 
It is easy to see that this system is satisfied if 
qi = v, ^ = 2a/^^(^ 
or 92 = 12 , ^ = 21 / 7 ^ 1 ( 91 ). 


(40) 


(41) 


Both solutions lead to the same physical solution due 
to symmetry of Eqs. m with respect to transposition 
of 9 i and 92 . For definiteness we shall take the second 
solution in m- Then the variable 91 oscillates in the 
interval < 91 < 1 ^ 2 , where 72 . 1 ( 91 ) > 0 and, hence, 
for /3 > 0, we have according to (l26l) two choices for the 
parameters u, z^i, 1 / 2 , 


P <vi < V2^ -fd <v < P 
-fd <Vi <V2< P, P> P- 


and their substitution into o yields the flow velocities. 
These formulae represent the periodic nonlinear wave 
which can be called the density wave, since the densi¬ 
ties oscillate in phase and in the small amplitude limit 
this wave reduces to the sound wave, which describes os¬ 
cillations of the total density p = Pi -\- P 2 - 

Let us consider the soliton limit when 123 —^ V 2 {m —>■ 

1 ): 


Pi = 


P2 = 


v + p 
2P 
i^-P 
2p 


P + V2- 
P-V2 + 


V2 - Vi 


cosh (^j 22 - viSP 
V2 - Vl 

cosh^(VJ22 - yiC) 


(46) 


where P < vi < V 2 and —P < v < p. These parameters 
can be expressed in terms of the constant densities at 
l?l 00 , 

Pw = ^(p+P){p+V2), P20 = ^{v-P){P-V2). (47) 

Solving this system with respect to i >2 and v gives 
122 = 2 (pio + P 20 — \/(pio — P 20 — 2/3)^ -b 4 piop 2 o) ) 

i 2 = - ^pio + P 20 + 'J (pio — P 20 — 2/3)^ -b 4 piop 2 o) ■ 

(48) 


The parameter P can be also expressed in terms of 
PiOi P 20 and constant flow velocities uio, U 20 at |^| —>■ 00 . 
From © and ([5]) we get 


0!i = Pio(uio — V), a 2 = P 2 o{u 20 — y) (49) 


and 


Pi = 2 (“10 “ T’lo + P 20 , 
P2 = 2 (^20 — yY + pio + P20, 


( 42 ) 


( 50 ) 
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FIG. 2: Dependence of the inverse width k of a dark-dark 
soliton on its velocity V for different values of the relative 
velocity between condensates. In all plots pio = P 20 = 0.5, 
uio = 0 and V is measured with respect to quiescent conden¬ 
sate: (a) U20 ~ 0, this is the case of well-known dark-dark 
solitons propagating through two still condensates; (b) plot 
corresponding to U20 = 1-5 illustrates a non-monotonic de¬ 
pendence; (c) at U20 = 2.0 the region of possible velocities 
splits with formation of two disconnected regions; (d) plot for 
U20 = 2.2 illustrates appearance of two disconnected regions 
of possible soliton velocities V. 


hence 


/3 = i [(«io - Vf - {u 2 o - F)"] . (51) 

To determine the last unknown parameter vi , we remark 
that Viete formula for the polynomial (P51) in our case 


gives Pi + 132 = vi + 2 (z /2 + v) and, consequently, 

The solution (l46l) exists if V 2 > I'l and this condition 
gives restrictions for the soliton velocity, 


Pw 


P20 


(V-Uio)'^ {V-U 2 o)‘^ 


> 1 . 


(53) 


Note that for the second choice in (1421) the condition 
1^2 > cannot be fulfilled. In the limiting case of 
one-component quiescent condensate (p 20 = 0, uio = 0) 
the condition (1551) reduces to the well-known fact that 
the soliton velocity is smaller than the sound velocity, 
V < Cs = i/pTo- 

Dependence of the inverse width k = ^V 2 — vi of soli¬ 
ton on its velocity V is given by 


■^|pio + P 20 — (^^10 — — {U 20 — 


+ \/[pw - P 20 - (wio - y)^ + (W 20 - + 4pio/02o} 


1/2 
(54) 


This expression can be also obtained by linearization of 
equations dH) with respect to small deviations around 
asymptotic densities {pk = Pko + Pfc) and seeking the 
solution of the linearized equations in the form p'^ oc 
exp(— k|^|). This calculation shows that the Appelrot- 
Delone class of solutions yields in the corresponding limit 
all soliton solutions with exponentially decaying tails 
around non-zero background densities pko ^ 0. 

Dependence (l54l) is illustrated in Fig. 2 for different 
values of the relative velocity of the BEG components. 
The remarkable new feature is that this dependence can 
be non-monotonic and for large enough values of the rel¬ 
ative velocity the region of possible values of the velocity 
y splits into two separated regions in sharp contrast with 
the one-component situation. The appearance of two re¬ 
gions of velocity can be illustrated graphically in the fol¬ 
lowing way. We introduce for convenience the variables 
X = V — uio, V = y — U 20 ; then the boundary of the 
region (ESI) is given by the equation 

- pioV^ - p2oX^ = 0. (55) 


Its plot is shown in Fig. 3 by a solid line and the ad¬ 
missible values of V are located inside this line (that 
is in the area including the origin of the coordinate 
system). If we fix the value of the relative velocity 
Uo = U 20 — uio = X — V, then the possible values of F 
correspond to points of the straight line located between 
its intersections with the curve (I55|) . Consequently, the 
splitting of the region of possible values of V correspond 
to such Uq that the straight line X — Y = Uq touches the 
curve ()55|) at the point where dY/dX = 1 or 


^y2 ^2y _ - P 20 X = 0 


(56) 
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(see Fig. 3). The system (|55]) and (l56t can be easily 
solved to give 


X = V-u,o = + 

Y = V-U2o = TP2i^\JpYo' + pli\ 


(57) 


and hence the critical value of the relative velocity is 
given by 


Uo = ip\i" + pli"f^^- (58) 



FIG. 3: Plot of the curve (1641) (solid line) and of the straight 
lines X — Y = Uo (dashed lines): (b) corresponds to a single 
region of possible values of V (see Fig. 2b); (c) correspond to 
the relative velocity at which the region splits into two regions 
(see Fig. 2c); (d) corresponds to two separated regions (see 
Fig. 2d). 


As in the preceding subsection, we see that the second 
equation is satisfied identically hy q 2 = P and the first 
equation dqi/d^ = — 2(/3 — (?i)-\/<?i — vi can be easily in¬ 
tegrated to give qi = j3 — [jd — vi ) cosh“^ WP — vi^). As 
a result we obtain the densities 


Pi = (i/-h/3) 


{P-vi)l{2P) \ 

cosh^(v'/3- 


P2 = {v 


{P-ui)/{2P) 

COsh^(-v/;9 — Vi ^) 


(61) 


which obviously correspond to the dark-bright soliton: 
the density pi has a dip at ^ = 0 and approaches to the 
background density po = v + /3 as ^ oo whereas p 2 
has a hump at ^ = 0 and vanishes as |^| —>■ oo. Let 
us relate the parameters of formulae (IFTI) with standard 
physical parameters for the soliton solution. To this end 
we define the inverse half-width k of the soliton by the 
equation k = i/(d — vi and introduce the ratio of the 
components densities at the center of the soliton 7 = 
(po — Pi( 0 ))/p 2 ( 0 ) = {P + P)/{P — P). Besides that we 
assume that there is no flow of the first component at 
infinity: Mio = 0. Then from Eqs. © and dS]) we find 
ai = -Vpo, 02 = 0, /3i = Po + /32 = Po - k^/2 and 

hence 


T 


7 “ 1 

y “ 27 ’ 


Vl = p - V = Pq 


7 -f 1 
(62) 

The dependence of the soliton’s inverse width on its ve¬ 
locity is given by the formula 


niy) 



(63) 


These formulae are equivalent to those found in Ref. 


M- 


C. Dark-bright soliton solution 


D. Legendre-Jacobi class of solutions 


If one of the background densities vanishes (say, p 2 o = 
0 ), then the so-called dark-bright soliton solutions of the 
Manakov system are obtained (see, e.g., EEl). Here 
we show that this type of solutions is a specialization 
of general solutions of the Kowalevski equations when 
the polynomial TZ{q) has two double zeroes. In this case 
the condition (l29ll is fulfilled if one of the double zeroes 
coincides with /3. Thus, we assume that 


— P dd qi < P and P dk q 2 dd V = v/i = v^. (59) 
Then the system (1^ reduces to 


2{P - qi){v 
2(92 - P){P 


91) V9i - = (91 - 92) 

92) ^92 - = -(91 - 92) 


dqi 

dC 

dq2 

dy 


(60) 


The general one-phase traveling waves described by the 
Manakov system can be illustrated by an easy numeri¬ 
cal solutions of the Kowalevski equations (1^51) . On the 
other hand, the analytical solution (l?7l) can be expressed 
in terms of Riemann 0-functions by the methods used 
already by S. Kowalevski (see and developed 

further in the algebr aic-g eometric approach to integrable 
equations (see, e.g., 0). This method was applied to 
the one-phase solutions of the focusing Manakov system 
in Refs, [l^, [l^. However, the resulting expressions are 
quite inconvenient for practical use. Therefore we shall 
confine ourselves here to a particular case, when the so¬ 
lution can be reduced to the much better known special 
functions (elliptic integrals) which permit one to under¬ 
stand the characteristic features of the solution in a much 
simpler way. Here we shall consider such a situation first 
noticed by Legendre and generalized by Jacobi [^ . 



























Let the zeroes of the polynomial TZ{q) be given by 
1/1 =c- 1 / 6 , i /2 = c- 1 /a, 1/3 =c, 

(64) 

i/4 = c+l, 1/5 = 0 + 1 / 06 , 

where 0 < 6 < a < 1 and the parameter c satisfies the 
conditions 

max{/3, 1/6 —/?} < c </? + 1/a. (65) 

so that qi and 52 oscillate within the intervals 

— P < Vl < qi < V2 < P, P < < q2 < Vi. (66) 


Let us assume for definiteness that at ^ = 0 we have 
qi{0) = c — 1 /a and q 2 = c (other choices of the initial 
conditions can be considered in a similar way). Then, 
introducing the variables 


■ 21,2 = qi .2 - c. 


(67) 


we represent the solution (I27p in the form 




zdz 

v/^ 


= ± 2 ^, 


where 


( 68 ) 


TZ{z) = z{l — z){l — abz){{l + az){l + bz)/{ab)‘^. (69) 


As Jacobi showed [ 2 ^, the integrals here can be calcu¬ 
lated in terms of incomplete elliptic integrals of the first 
kind. Since Jacobi did not provide the details of his 
method, this calculation is discussed briefly in Appendix. 
As a result, we obtain a particular solution of Eqs. (El 
in the form 


ki) + k2) - F{(P2, ki) - F{(P2, ^ 2 ) = 0, 

F{(pia, ki) - F{ipib, k2) - F{ip2,ki) + F{ip2,k2) 


= ±2 


(1 + a)(l + 6 ) 


ab 




(70) 


These equations determine implicitly the dependence of 
Zi and Z 2 , and, hence, of qi and q 2 , on in the interval 
of ^ until the first turning point is met (zi = — 1/6 or 
Z 2 = 1). After that the sign before the corresponding 
square root in the Kowalevski equations (l2^ must be 
changed and the replacement in the solution (l68ll 

n dz p dz 

\Jn{z) ^Jn{z) ^n{z) 

or 

dz dz _ n dz 

must be done with similar changes in the expressions 
dZQl). Making such changes at every successive turning 
point, we find the solution in any necessary interval of 
Substitution of resulting qi = zi + c and 52 = -22 + c into 
Eqs. (ITSl) yields the dependence of densities pi and p 2 on 
y Typical resulting plots are shown in Eig. 4. 





where 


+la = 


arcsm 


(l-|-azi)(l+bzi) 

zi(v/a—+ 6)2 


7r-arcsmW ^~'7'+^~+A\ —^ < 21 < 


(l-|-azi)(l-|-bzi) 
zi(+a—+6)^ 


1 ^ ^ 1 . 
’ +76 - - a ’ 


(71) 


/(I + azi)(l + 6zi) 1 1 

(/?ib = arcsm y-—- “t < 21 <-; 

V zi ya + V 6)2 6 a 


(72) 


P2 = arcsm i 


(1 + o) (1 + 6 )z 2 
(1 + az 2 )(l + bz2) 


, 0 < Z 2 < 1. (73) 


FIG. 4: (Color online) Plots for the solution (1701) of the 
Kowalevski equations for the values of the parameters a — 0.8, 
b = 0.4, c = 2 which correspond to vi = —0.5, 1/2 = 0.75, 
i /3 = 2.0, i /3 = 3.0, 1/5 = 5.125. The initial conditions are 
given by 51 ( 0 ) = A 2 , 92 ( 0 ) = A 3 , (a) Plots of qi and 52 for 
the interval of ^ corresponding to full cycle of q 2 -variable; 
dashed lines indicate the values of Vi, i = 1,2,3, 4; (b) plots 
of the components densities pi and p 2 and of the total density 
p = pi + P 2 . 


As one can see in Fig. 4, in the general solution the 
periodicity of the wave in space and time is lost and the 
wave pattern demonstrates quite complicated behavior 
as a function of 
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IV. CONCLUSION 

In this paper we have found the one-phase traveling 
wave solution of the Manakov system which describes 
evolution of two-component BEC. It is shown that in 
this case the Manakov system reduces to the equations 
which S. Kowalevski derived in her study of rotation of a 
heavy top in the discovered by her completely integrable 
case. We show that the previously found solutions of 
the Manakov system appear in this scheme as particu¬ 
lar cases. Besides that, new solutions are found which 
were either missed in previous analysis or cannot be ob¬ 
tained by more elementary methods when parameters of 
the solutions are chosen in such a way that the evolu¬ 
tion equations are greatly simplified. In particular, we 
have found a new dark-dark soliton solution for a two- 
component BEC with the non-zero relative motion of the 
components. This solution has very unusual dependence 
of the inverse width on the soliton’s velocity. In princi¬ 
ple, this can lead to new forms of dispersive shock waves 
evolved from initial step-like distributions of the compo¬ 
nents densities or velocities. 

For applications of the developed theory to the de¬ 
scription of the jMlarization wave patterns observed in 
the experiments the Whitham modulation theory 

for these waves has to be developed. Some par¬ 
ticular situations have already been studied in Refs. 
(genus-zero case) and [l^ (genus-one case for the limit 
P = 0). The results of the present paper demonstrate 
that the general one-phase solution is described by the 
fifth-degree polynomial TZ{q) whose zeroes as well as the 
wave velocity must be related in framework of the finite- 
gap integration method with the modulation parameters 
appearing in the Whitham theory of modulations of non¬ 
linear waves. Thus, the results obtained here provide the 
necessary step to development of the modulation theory 
which can be applied to description of dispersive polar¬ 
ization shock waves observed experimentally. Derivation 
of the Whitham equations is a difficult problem far be¬ 
yond the scope of this paper and we hope to consider it 
elsewhere. 
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Appendix 


The integrals Ii and are calculated with the use of the 
substitution 


1 -I- abz'^ = uz 


(A.2) 


or 



where u is a new integration variable. It is easy to see 
that the function (IA.3I) with the lower sign maps the 
interval I-|-a&<M<ooonO< 2 ;<l and with the 
upper sign maps the same interval on l/{ab) < z < oo. 
Then substitution (jA.3|) with the lower sign into Ii gives 
after simple manipulations 



(A.4) 


where 


u{z) = {1 + abz‘^)/z. (A.5) 

Elliptic integrals in (jA.4l) are transformed to the standard 
form by the substitution 


( 1 -I-a)(l-I-6) 

u = -- a — b. 


sin^ (fi 


(A.6) 


As a result we obtain 
ab 


h = 


'\/(l + ®)(1 + 


{F{ip,ki) + F{Lp,k2)} , (A.7) 


where F{if, k) denotes the elliptic integral of the first 
kind. 


ki = 


^/a- Vb 




k<2 = 


-/a -I- y/b 


(k.S) 


and is related with the upper limit of integration z by 
the formula 


. 2 

sm ip = 


(I -f o)(l -f b)z 
(1 -I- az){l + bz) 


(A.9) 


The integral is calculated by the same method and the 
result reads 


We have to calculate the integrals 

dz ,, zdz 


h = 




n{z) 


n{z) 


h 




(A.l) 


= + (A.io) 

The integrals I 2 and in (lA.IIl can be calculated with 
the use of the substitution 


u{z) = —(I -I- abz‘^)lz 


(A.ll) 
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or 


where 



+ 2'/ab ± 



which map the interval 2'/ab < u < a + b on the intervals 
— 1/5 < z < —Xj'/oh and —Ij'/ah < z < —1/a, corre¬ 
spondingly for upper and lower signs. This transforms I 2 
to 


^ ^ _ dM _ 

^ \J[u — a — b){u + 2\fal)){u -I- 1 -I- ab) 

^ ^ _dn_ 

^ \J [u — a — b){u — 2\/ah){u + 1 + ab) 

(A.13) 

These integrals are reduced to standard form of elliptic 
integrals by substitutions 


TT — arcsm 


Va = 


and 


(l+az)(l+bz) 


— ■r < Z < - 

^ y/ah^ 


arcsm 


(l-eaz)(l-|-bz) 
z{y/a—\/b)‘^ ' 


< z<-- 


Vob 


(A.16) 


= arcsm, 


' {1 + az){l + bz) 
z{^/a + ^/b)‘^ 


Similar calculation yields 


1 1 

-T<z<—. 
b a 


(A.17) 


M = a-|-5 — (a-|-5± 2-\/o5) sin^ ip. 

As a result we obtain 
ab 


(A.14) 


l2=- 


\/(l-|-a)(l-|-5) 


{F{pa,ki) + F{pb,k2)} , 


(A.15) 


^ 

^ V (1+®)(i+ 


{i^((/9a,fcl) - ■F’(V56,fc2)} • (A.18) 
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